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Abstract. The migration of profile sub-peaks identified in time-monitored optical emission lines of Wolf-Rayet 
(WR) star spectra provides a direct diagnostic of the dynamics of their stellar winds via a measured AvLOs/At, a 
line-of-sight velocity change per unit time. Inferring the associated wind acceleration scale from such an apparent 
acceleration then relies on the adopted intrinsic velocity of the wind material at the origin of this variable pattern. 
Such a characterization of the Line Emission Region (LER) is in principle subject to inaccuracies arising from line 
optical depth effects and turbulence broadening. In this paper, we develop tools to quantify such effects and then 
apply these to reanalyze the LER properties of time-monitored WR stars. We find that most program lines can 
be fitted well with a pure optically thin formation mechanism, that the observed line-broadening is dominated 
by the finite velocity extent of the LER, and that the level of turbulence inferred through Line Profile Variability 
(Ipv) has only a minor broadening effect in the overall profile. Our new estimates of LER velocity centroids are 
systematically shifted outwards closer to terminal velocity compared to previous determinations, now suggesting 
WR-wind acceleration length scales /3R* of the order of 10 — 20 Rq , a factor of a few smaller than previously 
inferred. Based on radiation-hydrodynamics simulations of the line-driven-instability mechanism, we compute 
synthetic Ipv for Ciii5696A for WR 111. The results match well the measured observed migration of 20-30 ms -2 , 
equivalent to /3R, ~ 20 Rq. However, our model stellar radius of 19 Rq, typical of an O-type supergiant, is a 
factor 2-10 larger than generally expected for WR core radii. Such small radii leave inferred acceleration scales to 
be more extended than expected from dynamical models of line driving, which typically match a "beta" velocity 
law v(r) = Voo(l — Rt/r) 13 , with /3 ~ 1 — 2; but the severity of the discrepancy is substantially reduced compared 
to previous analyses. We conclude with a discussion of how using lines formed deeper in the wind would provide 
a stronger constraint on the key wind dynamics in the peak acceleration region, while also potentially providing 
a diagnostic on the radial variation of wind clumping, an issue that remains crucial for reliable determination of 
O-star mass loss rates. 

Key words, line: formation - radiative transfer - stars: atmospheres - stars: early type - stars: mass loss - stars: 
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1. Introduction 

The hot, massive, luminous stars of spectral type O, B, 
and WR are understood to have strong stellar winds 
driven by the line-scattering of the star's continuum radi- 
ation (see, e.g., Puis and Kudritzki 2000). For OB stars, 
dynamical models based on the line-driving formalism of 
Castor, Abbott, & Klein (1975; hereafter CAK) and its 
modern extensions (Pauldrach, Puis, & Kudritzki 1986; 
Friend & Abbott 1986) yield predictions for the mass 
loss rate M and terminal speed Voq that are generally in 
quite good agreement with observationally inferred val- 
ues (Puis et al. 1996; Crowther et al. 2002). Indeed, for 
these relatively optically thin winds, the formation of line 
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profiles throughout the flow acceleration region yields a 
constraint on the full wind velocity law, commonly char- 
acterized through the phcnomenological "beta-law" form 
v(r) — Voo(l — R^/r) 13 , where i?* is the (hydrostatic) sur- 
face value of the radius r, and the velocity exponent /3 is 
inferred from line profile fits. Careful time- monitoring of 
optical emission lines (e.g. He II 4686; Eversberg, Lepine, 
& Moffat 1998) show that these winds are not strictly 
steady, but likely consist of quite extensive small-scale 
structure. This leads to an intrinsic, low-level, Line Profile 
Variability (Ipv), characterized by moving subpeaks that 
provide an even more direct diagnostic of the wind acceler- 
ation. For OB stars, both the single exposure line profiles 
and such Ipv subpeak monitoring suggest a velocity-law 
exponent (3 ~ 1 — 1.2, which is roughly consistent with 
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the value (3 ss 0.8 — 1 predicted from dynamical models of 
line-driving. 

Over evolutionary timescales, the cumulative depletion 
of the hydrogen-rich stellar envelope of such OB stars is 
understood to lead to the Wolf-Rayet (WR) stage, char- 
acterized by reduced or no hydrogen, with spectra domi- 
nated by broad emission lines of He and CNO processed 
material. For such WR stars, the very large mass loss rates 
inferred from the strength of these wind-broadened emis- 
sion lines are difficult to reconcile with a detailed, dynam- 
ical model based on the paradigm of radiative driving via 
line-scattering. Because the higher density makes the in- 
ner wind optically thick even in the electron scattering 
continuum, there are now no direct diagnostics of the ini- 
tial acceleration from the hydrostatic stellar core, making 
it difficult to constrain dynamical models. Moreover, the 
overall shape of the emission profile turns out to be rela- 
tively insensitive to even the acceleration scale in the outer 
part of the wind (see Sects. 3.1-3.2). 

For such WR winds, the monitoring of moving sub- 
peaks from Ipv thus provides a crucial diagnostic of the 
outer wind acceleration (Robert 1992). The most exten- 
sive previous analysis of such Ipv is that of Lepine & Moffat 
(1999; hereafter LM). By combining the measured sub- 
peak acceleration with an estimate of the characteristic 
velocity and range of the Line Emission Region (LER), 
they inferred wind acceleration scales w 20 — 80 Rq . 
If applied to typical OB supergiant radii i?» w 20 — 45 Rq 
(Crowther et al. 2002; Hillier et al. 2003), such accelera- 
tion scales imply a velocity index (/3 ~ 1 — 2), which is 
nearly compatible with the value predicted by wind mod- 
els. But both evolutionary (Langer et al. 1994; Schaerer 
& Maeder 1992) and spectroscopic (Hillier & Miller 1999) 
models of WR stars suggest a much smaller hydrostatic 
core radius, i?* w 2 — 5 Rq , for which the LM results im- 
ply a very large ratio between acceleration scale to core 
radius, formally (3 w 4 — 40, with a typical value (3 w 10. 
Over an associated factor 10 range in radius the radiative 
flux declines by a factor of 1/100, making it difficult to 
understand how such extended acceleration scales could 
be maintained by a radiative driving mechanism. 

This difficulty motivates our present reexamination of 
the basic inference of wind acceleration scales from Ipv, 
with particular emphasis on the assumed velocity and 
width of the LER. We first provide a simple overview 
(Sect. 2) of how a measured acceleration can be related to 
an characteristic acceleration scale through an associated 
velocity range. We next (Sect. 3) review the sensitivity of 
the overall shape of emission profiles to the wind and LER 
parameters, showing thereby that the observed emission 
line broadenings are consistent with relatively large and 
extended velocities for the LER, without need to invoke 
large turbulent speeds. We then (Sect. 4) show that such 
larger and more extended LER velocities imply substan- 
tially lower acceleration scales (/3i2* « 10 — 20 Rq), which 
are now much closer to being compatible with dynamical 
models that account for the multi-line scattering expected 
in dense, WR winds. In Sect. 5 we examine how well Ipv 
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Fig. 1. The radial distance to reach from a given scaled 
velocity wi to half-way to the terminal speed, (1 + wi)/2, 
scaled by /3i2*, and plotted vs. wi for various velocity in- 
dices (3 — 1/2 (lowermost curve) and (3 =1, 2, 4, 6, 8, and 
10 (overplotted curves). 

and associated subpeak accelerations can be reproduced 
through dynamical simulations of the intrinsic instability 
of line-driving. Finally, Sect. 6 summarizes our results and 
provides an outlook for future work. 

2. Dynamics of wind acceleration 

2.1. Acceleration length scale 

To gain a basic understanding of the factors influencing an 
inferred wind acceleration scale, consider first the simple 
case of constant acceleration a between a starting velocity 
vi and final velocity v^. By simple dynamics, the associ- 
ated length is 



si 



2n 



(1) 



For application to a stellar wind, let us scale the starting 
velocity by the terminal value, wi = vi/voc, and the length 
by a solar radius, 



si 

Rq 



1 



Wt 



2h 



(2) 



where 5 = (iRq/v^ represents a convenient scaling for 
the acceleration. Applying now the acceleration magni- 
tudes inferred by LM from WR Ipv, a » 0.001 - 0.005, we 
immediately derive associated length scales 



si 

Rq 



= 100 - 500 (1 - wf) . 



(3) 



If we now associate the quantity wi defined above with the 
characteristic fractional velocity of the inner edge of the 
LER, we see that the typical LER speeds assumed by LM, 
wi w 0.6 imply very large acceleration scales, si/Rq rs 
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60 — 300. By contrast, if we instead take wi « 0.9 (see 
Sect. 3.2, Fig. 4), we find s;/i? Q w 20 - 100. 

An alternative way to characterize the overall wind 
acceleration is to make use of the phenomenological beta- 
velocity law 



w(r) = v(r)/voa = (1 - R*/r) 



(4) 



for which the scaled acceleration vs. velocity becomes 

a = R Q ww' = ^(1 - w^fw 2 - 1 /?. 
R* 



(5) 



However, since w < 1, note that for moderately large /3, 
W i/P =e(W }nw ttl + (l/j3)]nw, (6) 



which implies 



(w lnw)^ 



(7) 



This very useful property, which was first derived by 
Lepine (1998), implies that a plot of acceleration vs. speed 
just scales homologously with the parameter (see Fig. 
8 below). For an assumed scaled velocity of the LER, wi, 
this now gives a scaled acceleration length, 



PR* _ (wi lnw;) 2 

Rf) a 



(8) 



Comparison with Eqn. show that the beta-law char- 
acterization of the acceleration scale differs from that 
for constant acceleration by the replacement of numer- 
ator (1 — wf) with (wj lntc;) 2 . Applying now the LM 
values for the velocity of the LER, wi « 0.6, we find 
(3R*/Rq « 20 — 90. By contrast, taking wi fts 0.9 gives 
@R*/R Q « 2 - 9. 

Note that these are much smaller than the simple, lo- 
cal scale given by Eqn. J5J. The reason for the difference 
is illustrated in Fig. 1, which plots the distance to reach 
to half the terminal speed as a function of initial speed, 
scaled by the characteristic scale f3R* , and overplotted for 
various values of the velocity index (5. The order unity val- 
ues for small initial velocities shows that /3i?* is indeed a 
characteristic acceleration scale for the inner regions of a 
beta velocity law; but at large velocities, the local accel- 
eration scales are actually much larger. 

2.2. Simple overview of radiative driving 

For a steady-state, spherical wind with velocity v(r) as 
function of radius r, the net flow acceleration a = v(dv / dr) 
results from the outward driving against the inward pull 
of gravity, 



dv 

V— = g r ad 

dr 



GM 



(9) 



where G and M are the gravitation constant and stellar 
mass. For WR and other very luminous stars, the general 
belief is that the outward driving is from stellar radiation 
interacting with both line and continuum opacity of wind 



material. For a star of luminosity L, the radiative flux 
is F = L/Anr 2 , implying that overall the radiative ac- 
celeration tends to have the same inverse-radius-squared 
decline as gravity. Since gravity moreover sets the over- 
all acceleration scale, it is convenient to define the gravi- 
tationally scaled accelerations v! = r 2 v(dv/dr)/GM and 
Trad = r 2 g rad /GM. 

For driving associated with continuum opacity k c , 
g ra d = k c F/ c, and so 



r c 



inGMc 



(10) 



In the common case that the continuum is dominated by 
electron scattering opacity n e (w 0.34 cm 2 /g), this is a 
constant that just reduces the effective gravity by a factor 
1 — r, with r the classical Eddington parameter. 

By contrast, for line opacity, the Doppler-shift deshad- 
owing of optically thick lines gives line driving a depen- 
dence on the flow acceleration, parameterized within the 
CAK formalism as 



fCu" 



(11) 



where < a < 1 is the CAK power-index, / is a complex 
function of velocity, radius, and acceleration that accounts 
for the finite size of the stellar disk (see CAK Eqn. (50)), 
and C is a constant set by line-opacity, luminosity, and 
mass loss rate. 

With these scalings, the equation of motion (9) can 
thus be written as 



u' = fCu ,a - (i - r) . 



(12) 



For the idealized case of a point-source star, f — I, this 
equation has no explicit radial dependence, implying that 
each term is just a fixed constant. A trivial radial inte- 
gration of the scaled acceleration u' from a static lower 
boundary stellar radius R* thus yields the velocity law 



v(r) 



; (i-i?*A) 1/2 . 



(13) 



where the terminal speed — v esc vw is propor- 

tional to the effective surface escape speed, v esc = 
V /2GM(1 - r)/i?». The critical CAK solution with the 
maximal mass loss rate further requires 



1 



du' 



(14) 



which when combined with Eqn. (12) implies a criti- 
cal acceleration u' c — C c — a/(l — a), yielding Uqo = 
v esc ^a/{l - a). 

For more complete models that take account of the 
finite stellar disk (Friend & Abbott 1986; Pauldrach et al. 
1986), the reduction of the line- force near the star (where 
/ < 1) lowers both the initial acceleration and the mass 
loss rate. But as / increases outward - with / — + 1 at large 
radii, where the stellar radiation indeed approaches the 
radial streaming from a point-source - the lower mass loss 
rate allows a much stronger acceleration (u' « 6 — 9!), thus 
yielding a faster wind terminal speed, Voalv esc = 2.5 — 3. 
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The resulting velocity law is no longer analytic, but can 
still be well fit by a general 'beta-law' form Q , with now 
13 w 0.8 - 1 (Pauldrach et al. 1986). 

For WR stars, the importance of line overlap severely 
complicates the derivation of self-consistent dynamical so- 
lution, but it is possible to extend the basic CAK for- 
malism if one makes the simplifying assumption that the 
lines are uniformly distributed throughout the stellar spec- 
trum (Friend & Castor 1983). Using this approach, multi- 
line scattering models based on Monte Carlo simulation 
(Springmann 1994) and on non-isotropic diffusion anal- 
ysis (Gayley, Owocki, & Cranmer 1995) both show that 
effective trapping of radiation in the expanding envelope 
can lead to a more extended acceleration, with an effective 
velocity index up to f3 ~ 2 in the outer wind. 

An overall point from such dynamical models based 
on line-driving is that it is quite difficult to obtain a ve- 
locity law with large velocity index, i.e. j3 > 2. This is 
in part because line-driving itself depends on the velocity 
gradient, and so becomes inefficient in the extended, slow 
accelerations of a large (3. Another fundamental factor re- 
gards the inverse-radius-squared fall-off of the radiative 
flux, F ~ 1/r 2 . Once the mass loss rate is tuned to allow 
the net outward force to be strong enough to overcome 
gravity near the stellar surface, the similar radial scaling 
of the radial flux and gravity tends then to make the net 
outward acceleration scale with the local inward acceler- 
ation of gravity. As shown above for the point-star CAK 
model, this leads to a velocity law with (3 = 1/2, and 
the terminal speed scaling with the surface escape speed. 
Corrections due to finite-disk factor or multi-line scatter- 
ing can extend the acceleration somewhat, but only up to 
velocity indices of order (3 « 1 — 2. Overall then, the very 
large velocity indices or associated acceleration scales in- 
ferred from apparent acceleration of subpeaks in WR Ipv 
are quite difficult to explain from dynamical models; this 
provides a key motivation for the present reexamination 
of these inferred extended acceleration scales. 



3. Line formation 

3.1. Discrete LER without broadening or stellar 
occulta tion 

To provide a basis for the analysis of Ipv, and its use 
as a diagnostic of flow acceleration, let us first examine 
how a mean, time-steady, line-profile depends on vari- 
ous wind properties. We again assume a spherically sym- 
metric, monotonically expanding wind described by the 
standard [3 velocity law, expressed in dimensionless form 
w(r) = i>(r)/uoo = (1 — i?*/r) /3 (Eqn. 4). As a simple 
first case, in this subsection we ignore the complications 
of stellar occultation and line broadening (i.e. by thermal 
or turbulent motions), and take the LER to have sharp 
boundaries defined by the lower and upper velocities v; 
and y u . The subsequent subsections will examine the ef- 
fects of modifying these simplifications. 



Following Owocki & Cohen (2001, Eqn. 7), and assum- 
ing optically-thin emission 1 , the total emitted luminosity 
at scaled wavelength x = (X/Xq — tyc/v^ is given by, 

LtMn(x) = ^ 2 C J dpj^ drr 2 p 2 {r)8{x + pw{r)), (15) 

where C is a constant that depends on the atomic process 
of emission (i.e. presently recombination). Performing the 
integral over p, gives 

rru r 2 p 2 {r) 

Lthin(x) oc / drH[w{r) - \x\\ -—— , (16) 

Jn w{r) 

where H is the Heaviside function (unity for positive argu- 
ment, zero otherwise). Note that from a shell at radius r, 
the contribution to the profile falls within the wavelength 
limits ±u>(r). 

Eqn. 1|16|) applies in the case when both the line and 
continuum are optically thin. For the case of an optically 
thick line, the relative escape of photons along direction 
cosine p is weighted by 

1 + cr(r)p 2 

9M= i + *{r)/r (17) 

where <r(r) = dlnv/dhir — 1. Integration over angle then 
gives (through the Dirac delta function) a corresponding 
weighting to the local emissivity, yielding for the optically 
thick emitted luminosity 

Lthick(x)<x f U drH[w{r)-\x\] T -^g{\x\/w{r),r) .{!%) 

In the optically thin case Q16[l . the wavelength x en- 
ters only via the Heaviside function; thus the emission 
profile from each localized, narrow shell (extending from 
r to r + dr) is constant within the wavelength bounds 
| a; | < w{r) defined by the scaled velocity w(r) of that 
shell. By contrast, in the optically thick case 1)18(1 . there 
is an additional wavelength dependence via the relative 
escape function g; this gives a curvature to the shell emis- 
sion profile, concave up in the inner wind where cr(r) > 0, 
and concave down in the outer wind where cr(r) < 0. 

These properties of line formation are illustrated in 
Fig. 2. The bottom panels show the individual flux contri- 
butions when the LER is divided into seven discrete, finite 
shells spanning equal velocity bins. In the optically thin 
case (left column), the contributions are frequency inde- 
pendent, yielding flat-topped, box functions. Progressing 
from higher to lower velocity shells, these contributions 
narrow while their emission level increases, following the 
increased density in the slower, inner wind regions. The 
middle panels show a cumulative stack of these individual 
shell contributions, thus approximating the total emission 
profile of the upper panels, which are obtained by a full 
radial integration over the continuum of differential shells 
that span the LER. Note that the inner velocity of the 

1 Here and throughout our analysis, we ignore continuum 
absorption. 
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Fig. 2. For optically-thin (left) and optically-thick (right) lines, the bottom panels show the emission contribution 
from seven discrete subshells with equal velocity width spanning an LER between 40% and 80% of Voo (with (3=1). 
The middle panels then show a stacking of these subshell emissions to build a cumulative line profile that approximates 
the spatially integrated emission profile from the full LER, as given in the top panel. 



LER sets the width of the flat top, while the outer ve- 
locity sets the maximum width of the profile down to the 
continuum level. 

For an optically thick line (right column), the forma- 
tion follows analogously, except that now the individual 
shell emission profiles can exhibit a net curvature within 
the box limits. In this particular example, wherein the 
LER spans the region where a(r) goes from positive to 
negative, the individual shell contributions go from con- 
cave up to concave down; the net cancellation of the two 
opposite senses of concavity gives the final profile a nearly 
flat-top form that is remarkably similar to the optically 
thin case. More generally, for optically thick lines formed 
predominantly in cither the inner (a > 0) or outer (<r < 0) 
regions, the total resulting profiles can exhibit either a 
double-peaked (concave up) form, or a rounded (concave 
down) form. 

In the limiting case of a constant-speed outflow, we 
have a = — 1 and g(p,r) = 3(1 — /i 2 )/2, yielding then upon 
angle integration a scaling 1 — x 2 for the (/-factor in Eqn. 
1)18(1 . which thus results in a parabolic shaped emission 
profile (see Ignace & Gayley 2002). 

Figure 3 compares optically thin line profiles computed 
for models with the same velocity range for the LER, but 
with much different velocity exponents, namely (3 = 1 and 
(3 — 4. Although these two values represent markedly dif- 



ferent spatial variations of the velocity, with much slower 
overall acceleration for the (3 — 4 case, the resulting line 
profiles are remarkably similar. This illustrates that the 
mean profile shape is quite insensitive to the wind acceler- 
ation, implying that inversion-based methods for inferring 
the wind velocity law from the profile shape (Brown et al. 
1997) are nearly degenerate with respect to the length 
scale of the velocity gain. 

3.2. Gaussian LER including stellar occultation 

Following NLTE model atmosphere results (Hillier 1988), 
as well as the convention from the previous Ipv analysis 
by LM, let us now specify the extent of the LER to follow 
a Gaussian form 

f(r) = exp(-[(«(r) - v e )/Av e } 2 ) , (19) 

where Av e is the velocity range for emission about the 
velocity centroid v e . For an outer wind velocity centroid 
v e = 0.9doo (and a relatively narrow velocity range of 
Av e = 0.1i>oo), Fig. 4 illustrates the resulting emission 
profiles for optically thin (dashed) and optically thick 
(solid) lines, comparing again a fast {(3 = 1; upper panel) 
and a slow {(3 = 4; lower panel) acceleration. The profiles 
are again all quite insensitive to the flow acceleration. Note 
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Fig. 3. Synthetic emission profiles for an optically-thin 
line forming between 30% and 90% of the terminal veloc- 
ity but for two different velocity laws: the broken (solid) 
line corresponds to a steep (shallow) acceleration with 
(3 — 1 ((3 = 4). Apart from a modified maximum line- 
profile flux (not known a-priori), the dependence of the 
line-wing shape on (3 is weak, i.e. ca. 10% flux-difference 
for an extended LER (shown here) and decreasing for nar- 
rower LER. Inversion techniques (Brown et al. 1997) yield 
a very uncertain information on the wind velocity law. 
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Fig. 4. Study of the impact of line optical depth effects 
and wind velocity law (/? = 1, 4) for a synthetic line whose 
LER is characterized by v e = .9«oc and Av e — .lwoc- 
Notice the rounding of the line-profile which shifts its 
HWHM velocity .3«oc below the adopted v e . Note also 
the negligible impact of [3 on the line-profile shape. 
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Fig. 5. Same as Fig. 4, but for an inner wind line with 
v e = .3voo and Av e = -l^oo- The double-peaking discussed 
in Sect. 3.1 is strong, as well as a profile asymmetry re- 
sulting from disk-occultation. This latter effect, observed 
in O supergiants, would permit the constraint of the wind 
acceleration just above the hydrostatic photosphere. 



that the Gaussian form of the LER now imparts a round- 
ing to the edges of the flat portion for the thin lines, but 
for thick lines the already curved form of the profile makes 
this additional rounding less noticeable. The vertical dot- 
ted lines mark the velocity centroid. Note that in both 
cases this lies outside the width at half-maximum, which 
is commonly used in observational analyses as a rough 
measure of the LER velocity centroid; for thick lines, such 
use of half-maximum leads to a substantial underestimate 
of the characteristic velocity of the LER. This can have 
important implications for the interpretation of Ipv, as we 
discuss below (Sect. 4). 

Figure 5 illustrates corresponding profiles for a line 
formed in the inner wind, with v e — 0.3woo. For the opti- 
cally thick lines, the double peaked profiles are the result 
of the concave-up form of the subshcll emission, as dis- 
cussed in Sect. 3.1; the marked Blue/Red asymmetry seen 
in the (3=1 case arises from occultation of the redshifted 
emission from the hemisphere behind the star. Occultation 
also has a smaller, but still noticeable effect on the red 
edge of the optically thin line, giving the overall profile 
a modest blueward shift. By contrast, occultation effects 
are much reduced in the (3 = 4 models, simply because 
the corresponding radius for the fixed velocity centroid 
v e = 0.3woo is now much further away from the star. 

Such B/R asymmetry is commonly observed in O su- 
pergiants, and so provides an useful constraint on the ac- 
celeration and velocity law (Puis et al. 1996; Crowther et 
al. 2002). 
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Fig. 6. Optically thin emission profiles for models with 
(3=1 and velocity centroid v e — 0.7voo, and three combi- 
nations of LER velocity width Av e and turbulent velocity 
^turb- For cases without turbulence («t U rb = 0), the solid 
and dashed curves show results for a wide (Av e = 0.3Vao) 
and narrow (Av e — O.lvoo) LER. The dot-dashed curve is 
for a very narrow LER (Av e = 0.05woo) but with substan- 
tial turbulence, fturb = ^-^ v oo- The vertical dotted line 
marks the velocity centroid v e . See Sect. 3.3 for further 
discussion. 

3.3. Effect of emission profile broadening 

To account for the effects of finite broadening (by ther- 
mal or micro-turbulent motions), we next apply the ray- 
integration method of Castor (1970; based on the Sobolev 
approximation (Sobolev I960)), with the narrow (delta 
function) emission now replaced by a Gaussian of velocity 
width wturb, 

<j>(x) <~ CXp [-(OToo/wturb) 2 ] ■ (20) 

Figure 6 shows the optically thin emission profiles for 
models with a common velocity centroid (v e = 0.7woo), 
but with three distinct combinations of LER width and 
turbulent velocity The solid and dashed curves compare 
a narrow (Av e = O.lwoo) and wide (Av e = 0.3foo) LER, 
without any turbulence (vturb = 0). Note that the lat- 
ter gives a centrally peaked profile with very extended 
line wings; again the half-maximum width thus gives a 
severe underestimate of the LER velocity centroid. The 
dot-dashed curve shows the profile of a very narrow LER 
(Av e = 0.05woo), but now with a substantial turbulent ve- 
locity («turb = 0.2^00); this turbulence further smoothes 
the edge of the profile, but does not substantially increase 
the overall profile width, even though such a turbulent ve- 
locity is at the upper limit of what is inferred in observa- 
tional analyses (Robert 1992; LM). This latter point plays 
an important role in our revision of the LM analysis of Ipv 
(see Sect. 4). 
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Table 1. Parameters used in profile fitting. The modeled 
line is Hcii541lA for WR 136 and WR 138, and Cm5696A 
for all other cases. To substantiate the terminal veloci- 
ties used in this study (first column), we compare values 
from LM (taken originally from from Prinja, Barlow & 
Howarth (1990), labeled a), and from various other stud- 
ies (b: Ecnncns & Williams 1994; c: Hillier & Miller 1998; 
d: Crowthcr & Smith 1996; e: Howarth & Schmutz 1992). 
Note that all values typically agree within ca. 100 km s _1 . 
The last column denotes whether the fit requires an opti- 
cally thick or thin line. 
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We can conclude from this section that a detailed nu- 
merical approach is required for the determination of the 
LER characteristics of broad and rounded profiles, since 
these profiles may be subject to significant optical depth 
effects. Observed WR emission profiles showing a central 
peaking and broad wings reflect a line formation over an 
extended LER. Because they do not show any double- 
peaking, these lines are most likely optically-thin. The ve- 
locity at half-maximum is a poor indicator of the LER 
velocity centroid, due to the bias of profile flux in favor 
of inner wind regions. Broad flat top profiles are an un- 
ambiguous sign of optically thin line formation close to 
terminal velocity (otherwise they would be narrow with 
minute flat top portion) and over restricted ranges in ve- 
locity (otherwise they would be peaked) . These latter sig- 
natures are the preferred diagnostics for Ipv studies. 

3.4. Application to line profiles observed at a time 
snapshot 

We are now in a position to use the model described above 
to assess the LER characteristics of time-snapshot obser- 
vations of the recombination lines Heii541lA (WNs) and 
CIII5696A (WCs) for seven WR stars: WR103, WR111, 
WR135, WR136, WR137, WR138 and WR140. This 
sample was chosen to match that used by Robert (1992) 
for her study of Ipv in WR stars, but the short-time ex- 
posures (ca. 30 min.) we use here were obtained indepen- 
dently by P. Crowthcr (p. a). (We exclude WR40 from 
our sample because its Hen line shows a narrow P-Cygni 
profile for which our method does not apply.) Let us reca- 
pitulate the model assumptions. We adopt a spherically- 
symmetric monotonically-expanding outflow. Line emis- 
sion occurs via recombination, thus possesses a density- 



8 Luc Dessart and S.P. Owocki: Inferring hot star wind acceleration from Ipv 



WR103 CIII5696 WR1 1 1 CIII5696 WR1 38 Hell541 1 




Fig. 7. Observed profiles (solid curves) and best fits (dot-dashed curves) in a chosen emission line (Hen 5411 A or 
Cm 5696A) for six WR stars discussed in LM. The intensity has been rescaled so that the continuum is at zero and the 
maximum intensity is at unity. The wavelength is in velocity units, scaled to the terminal velocity of the corresponding 
line, as given in Table 1. 



square dependence (see Eqn. 15). No continuum opacity 
is included but allowance is made for variable line-opacity. 
Turbulence is treated as in Eqn. (20). In Hillicr (1988), a 
detailed description of the ionization stratification is made 
for WR stars in general, showing that it is indeed realistic 
to use such a height segregation for recombination lines in 
WR atmospheres (see also Dessart et al. 2000). There is 
no explicit dependence of the emissivity on temperature, 
but it implicitly enters our treatment via this height seg- 
regation of wind ionization (Eqn. 19). We show our fits 
to single-exposure line profiles for these objects in Fig. 7, 
which shows the observed profiles (solid curves), together 
with our model fits (dot-dashed curves). (For clarity, we 
omit WR 135 from the figure; its profile resembles that of 
WR111). The model parameters are presented in Table 
1. The best fits are obtained by "eye-ball" estimate after 
iteration from a first guess, modifying the various parame- 



ters: line opacity, centroid and extent of the LER, without 
need to include any turbulence. 

Note first that only for one case - WR 136 - is it nec- 
essary to invoke line optical depth effects; this contrasts 
with the claim that such effects can be responsible for 
outstanding fit discrepancies (LM). For candidates other 
than WR136, we can identify some lines with extended 
line wings - for which we infer a broad LER (WR103 
and WR 138) - and others with near vertical line wings 
- for which we infer a narrow LER (the rest of the sam- 
ple of stars); overall, the velocity width spans the range 
Av e = 0.1 - 0.3^00. 

But a particularly important result is that the veloc- 
ity centroids of the LERs are all substantially higher than 
inferred by LM, in the range v e = 0.65 — 0.9woc. Indeed, 
when these centroids are coupled with the inferred widths, 
i.e. v e + Av e , we see that the line formation region in all 
case extends essentially all the way to the terminal speed, 
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. This fact is manifest by the close correspondence be- 
tween the maximum base width of these optical emission 
lines and the terminal speeds inferred by UV resonance 
lines. 

Note that our detailed modeling is capable of fitting 
line profiles with a very high level of accuracy since minor 
changes to model parameters lead to very different syn- 
thetic line profiles; the fit quality is markedly degraded 
if either v e or Av e are changed by even a few percents. 
The central part of the synthetic profiles for WR 103 and 
WR138 show minute variations due to numerical inac- 
curacies, resulting from the strong emissivity gradient at 
the base of the LER, while that of WR 136 shows a slight 
Blue/Red asymmetry due to disk occultation effects (Sect. 
3.2 and Fig. 5); In the present situations, such low- level 
irregularities, either real or from numerical error, would 
disappear with even a modest turbulence. 

With such information on the velocity in the formation 
region of emission lines with observed Ipv, we can now 
move on to using the wind acceleration derived from the 
migration of profile sub-peaks to infer the overall scale of 
the wind velocity law, as characterized by the value of 
(3R*. 

4. Application to inferred acceleration 

For these program stars, LM conducted an extensive anal- 
ysis of the Ipv observed by Robert (1992), inferring wind 
acceleration from the migration of profile sub-peaks, based 
on a phenomenological Discrete Wind Emission Element 
(DWEE) model. Using the inversion technique of Brown 
et al. (1997), LM estimated the LERs to extend over a nar- 
row velocity range (Av e « 0.05t>oo) located relatively deep 
in the wind (v e = 0.4 — 0.8uoo), with turbulence claimed 
to be the main profile broadening mechanism. 

Such narrow LERs warranted taking a constant wind 
acceleration over its limited range, a key assumption for 
the degradation function method used to infer wind ac- 
celeration magnitudes. Combining these with the inferred 
LER velocities, LM concluded that the associated flow 
acceleration must extend over very large length scales, 
w 20 — 80 Rq. These are much larger than the 
estimated hydrostatic core radii of WR stars (< 5 Rq, 
Hillicr & Miller 1999; Dessart et al. 2000, Crowther et 
al. 2002), a result that challenges the standard radiative 
driving paradigm for these WR winds (Sects. 1-2). 

By contrast, the LER characteristics derived here dif- 
fer markedly from those inferred by LM, with systemati- 
cally higher velocity centroids and larger widths, generally 
extending the outer range of the LER to nearly the ter- 
minal speed. Adopting the same acceleration magnitudes 
derived by LM, this now implies much more modest length 
scales for the wind acceleration. This follows simply from 
the fact that line formation in regions close to terminal 
speed can be expected to exhibit relatively low accelera- 
tion magnitudes. 

These central results are illustrated in Fig. 8, wherein 
the left panel reproduces Fig. 14 of LM, while the right 



panel reconfigures the LER location and width according 
to the results given in our Table 1. Comparing the left to 
the right panel, note that the acceleration values are the 
same, but the LER velocity centroids (crosses) are shifted 
rightward (by ca. 0.1— 0.3woo) and the velocity ranges (hor- 
izontal bars) are substantially broadened (typically more 
than doubled). As in LM, the overplotted curves in both 
panels show the associated variation of acceleration vs. 
speed for models with various acceleration scales , la- 
beled in units of Rq. 

Note in particular that our LER parameters in the 
right panel are now consistent with substantially smaller 
acceleration lengths (10 — 20 Rq) than implied by the LM 
values in the left panel (20 — 60 Rq). These are still quite 
large, but perhaps now within reach of what radiative 
driving can achieve if one accounts for multi-line scatter- 
ing (Springmann 1994; Gayley, Owocki, & Cranmer 1995). 

The origin of this difference lies solely in the new pa- 
rameters derived for the LER velocity location and width. 
LM cited turbulence as the dominant source of broadening 
of WR line profiles; but it seems that with LM's adopted 
narrow and deep LER, line profiles with extended wings 
can only be fitted by invoking very large turbulent veloci- 
ties, indeed much higher than those determined from their 
own Ipv analysis. 

Figure 9 illustrates these points for the sample case 
of the Heii541lA line from WR 138, comparing the ob- 
served profile (solid curves) with various models (dot- 
dashed curves). The top panels show our models with 
no turbulence (left) and with the Ipv-derrved turbulence 
(right); both fits are quite good, though the turbulence 
introduces a smoothing and a moderately improved agree- 
ment. The bottom panels show LM models with this same 
Ipv turbulence (left) and with the much larger turbulence 
needed to fit the observed profile (right); the former is the 
turbulence value cited by LM, but we clearly see that this 
does not provide an acceptable profile fit when combined 
with their inferred LER characteristics. Instead, a quite 
unrealistic turbulence is needed to fit the line. 

In the present models with the line emission formed 
over a range of velocities, the emission profile is rounded 
even without turbulence (Fig. 9, panel A), and so the 
addition of a modest turbulence leads to only a modest 
additional broadening without much affecting the overall 
emission profile shape (cf. panels A and B of Fig. 9). By 
contrast, in the LM models, in which the LER is assumed 
to be at nearly a single velocity, the intrinsic profile has a 
more flat-topped, box shape, even with the modest level 
of turbulence (i.e. Wturb = 220 km s _1 ; Fig. 9, panel C). 
Because of the relatively steep gradient in emission in go- 
ing from the core to wing, the stronger level of turbulence 
(«turb = 400 km s _1 now has the general effect of redis- 
tributing flux from near the edges of the flat-top into the 
line wings; the overall effect then is to broaden the base 
(wing) of the profile, but round-off and narrow the emis- 
sion core. An essential point is that the narrow range of 
velocity assumed by LM for the LER requires just such a 
large and physically dubious level of turbulence to produce 
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Fig. 8. (Left) A reproduction of Fig. 14 of LM. This figure depicts the local wind acceleration a e (determined from a 
degradation function analysis) at a given velocity v e in the stellar atmosphere for seven stars analyzed in LM (crosses) . 
The vertical error bar reflects an inaccuracy in the determination of the acceleration from the DWEE model of LM, 
while the horizontal error bar reflects the extent of, i.e. not an uncertainty in the position of, the LER in velocity 
space. The over-plotted curves are the theoretical a e versus v e for different values of /3i?*. (Right) Id. but with the 
LER characteristics determined in this work, i.e. crosses are shifted rightward by .05 — .25voo compared to LM, with 
considerably broader velocity ranges of the emission sites. 



the observed width and form of the line profile, whereas 
this can be produced naturally without invoking turbu- 
lence if one allows for the LER to span a modest range in 
outward flow speed. 

This exercise supports the view (Hillier 1988) that WR 
line broadening results primarily from the finite extent of 
the LER in velocity space, and not from turbulent broad- 
ening. These changes relative to the LER properties used 
by LM may not have much effect on the acceleration mag- 
nitudes derived in their Ipv analysis (see Sect. 5), but they 
do imply a substantial downward revision (by a factor of 
two or more) in their quoted acceleration scales, i.e. from 
20-80 R Q to 10-20 R Q . 

5. Wind acceleration inferred from synthetic Ipv 

The migrating sub-peaks observed in Ipv are understood 
to stem from the variable emissivity of the structured 
wind accelerating through the LER (Sect. 1). Lucy (1984) 
and Owocki & Rybicki (1984) proposed that this wind 



structure could result from the strong instability to small- 
scale perturbations of the line-driving mechanism. The 
non-linear evolution of a line-driven wind subject to such 
an instability was studied for the first time in a ID 
radiation-hydrodynamics computation by Owocki, Castor, 
& Rybicki (1988), showing that the outflow relaxes into 
a highly-developed non-monotonic radial velocity field, 
whose associated shocks lead to radially-confined slow 
high density compressions separated by extended fast rar- 
efied regions. We present these properties for one such 
simulation of an O-star wind (with suitable parameters for 
the supergiant star ( Puppis) in the top panels of Fig. 10 
for the radial velocity (left) and density (right) as a func- 
tion of radius. When plotted in a Lagrangian coordinate 
(bottom panels, Dessart & Owocki 2002b), the bulk of the 
material follows a velocity, which, when time-averaged, 
matches to within a few percents the CAK velocity law 
with a /3-exponent value of 0.8. In other words, only a 
negligible fraction of the total wind mass departs from 
such a CAK velocity, contributing relatively little flux to 
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Fig. 9. Top: Observed Hcii541lA profile of WR138 
(solid curves), compared with models using the parame- 
ters of the present study (dot-dashed curves) , without (A) 
and with (B, uturb = 220 km s _1 ) turbulence. Bottom: 
Same but with the parameters of LM (C), and with a tur- 
bulent velocity twice the value determined from the Ipv 
analysis (D, w tU rb = 400 km s _1 ). 



lack of emission profile variations in the line center re- 
gion. Moreover, Dessart & Owocki (2002b) interpreted the 
velocity-scale of variable profile sub-peaks at line-center 
as arising from the lateral geometrical broadening of wind 
structures. They showed that with a patch-size of 3 deg, 
one obtains a satisfactory fit to the observed sub-peak ve- 
locity width at line center as well as the fluctuation mag- 
nitude of observed Ipv. Thus, the same value of 3deg is 
used here. 

In the following two sections, we use the Patch method 
together with the hydrodynamical simulations of the ra- 
diative instability presented in Fig. 10 to compute syn- 
thetic Ipv of WR and O stars. Strictly speaking, such hy- 
drodynamical simulations only apply to winds driven via 
single line-scattering, i.e. O star winds, and not to WR 
winds whose larger mass loss rates require the treatment 
of multiple line-scattering (Gayley, Owocki, & Cranmer 
1995). Nonetheless, both are understood to be line-driven 
and believed to share similar dynamical properties, dif- 
ferences arising from the height location of line-emission: 
all O-star optical diagnostics originate close to the hy- 
drostatic core radius (below ca. 2-3 R», Hillicr ct al. 2002) 
while those of WRs subsist out to large radii. It is via such 
a height segregation that we now study such a synthetic 
Ipv for WR and O star emission lines. 



5.1. Measured sub-peak acceleration from synthetic 
Ipv: WR-star case 



(density-squared dependent) emission profiles. Thus, this 
property justifies the use of a standard /3-velocity law for 
the wind structures, as done in LM for the modeling of 
Ipv. 

However, because the lines used by LM have a finite 
rather than a narrow LER (Sect. 4), assuming a fixed 
wind acceleration may lead to systematic errors in the 
derived acceleration magnitudes. Indeed, in the case, e.g., 
of a wide LER, most of the profile flux originates from 
the low- velocity high-density regions (Sect. 3), which for 
a fully clumped wind, tends to favor the variability signal 
of those deep wind regions. In this section, to circumvent 
such a kind of biases, we present a more direct compar- 
ison between observed and synthetic Ipv. To do so, we 
compute the optically-thin line emissivity arising from a 
3D wind assumed to be composed of independent star- 
centered cones characterized by a fixed angular extent. 
The crux of this so-called "Patch" method (described in 
great details in Dessart & Owocki (2002a), hereafter DO) 
is to adopt, for each cone, a purely ID wind layout as com- 
puted by radiation hydrodynamics simulations of the ra- 
diative instability. A representative snapshot of the height 
variation of the velocity and density is shown in Fig. 10, 
both in Eulerian and Lagrangian coordinates. Thus, while 
the method retains some ID features as far as the gas 
is concerned, it allows us to introduce a complete wind 
anisotropy, which, as illustrated in DO, is required by the 



We stressed in the previous sections the importance of us- 
ing a reliable LER characterization for the interpretation 
of Ipv. In particular, we found that in general, the lines 
studied by LM formed over a finite LER that extended 
out to terminal velocity. Hence, here, for all profile com- 
putations, we set the outer LER radius at the maximum 
radius of the hydrodynamical simulation of 20-R*, corre- 
sponding to the terminal velocity of the flow within a few 
percents. We then vary the inner radius of the LER to 
reproduce the range of observed widths of Ipv diagnostics, 
covering the values n — 2.5,5,10 and 15i?». We show a 
sequence of synthetic grayscale time-series in Fig. 11, with 
increasing 77 value from top to bottom, and left to right. 

In the top-left panel, we see the presence of sub-peaks 
with a range of acceleration magnitudes. Around line of 
sight velocities of 1000 km s _1 and -1000 km s _1 , we can 
estimate values of ca. 50 km s _1 /hr, which due to a line 
of sight angle cosine weight of about a half, correspond to 
an higher intrinsic wind acceleration of ca. 30 ms~ 2 . In 
the line wing, the sub-peak migration magnitude is sys- 
tematically lower, ca. 20 km s _1 /hr and since these pro- 
file patterns arise from wind emitting structures advert- 
ing directly towards us, the sub-peak acceleration corre- 
sponds directly to the intrinsic wind acceleration, i.e. of 
ca. 8 ms~ 2 . This is consistent with the fact that the LER 
is wide, covering a large range of intrinsic wind acceler- 
ations. Despite additional geometrical effects, this leads 
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Fig. 10. Single time-snapshot from the hydrodynamics simulation of the radiative instability, showing 
(left column) and log-density (right column), in Eulerian (top) and Lagrangian (bottom) coordinates. 
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to a decreasing sub-peak migration magnitude from line 
center to line wing. 

As we move the inner edge of the LER to larger radii 
(subsequent three panels), the fast migration of sub-peaks 
in the flat part of the profile is at first less conspicuous 
(top- right panel) , before disappearing altogether for outer 
wind forming lines (bottom panels). In these latter cases, 
the LER covers a narrow range of velocities and is well 
characterized by a unique wind acceleration directly mea- 
surable from line-wing subpeaks. Contrary to the case of a 
deep forming line, the profile migration now increases from 
line center to line wing. Its magnitude is small, consistent 
with the fact that the LER now probes wind regions that 
are within 10% of the terminal velocity. 

Interestingly, the simulation corresponding to a line 
whose LER starts at 2.5 R* (top-left panel) repro- 
duces well both the temporal evolution as well as the 
time-averaged line emission observed for CIII5696A for 
WR111. In particular, the migration magnitude of ca. 30 
m s~ 2 is in good agreement with values given in LM as well 
as the observed sub-peak migration in the flat part of the 
profile. Note that the migrating subpeaks appearing in the 
line wing arc unobserved, a likely consequence of the fact 
that ID radiation hydrodynamics simulations underesti- 
mate the magnitude of radial velocity dispersion. But the 



more dramatic drawback from our model is that we use a 
stellar radius of 19 R*, while the hydrostatic core radius of 
WR 111 is expected to be a factor of 2-10 smaller (Hillicr 
& Miller 1999). Hence, with such small core radii, the 
inferred WR wind acceleration length scale is still signifi- 
cantly more extended than expected from radiation driven 
wind theory which at present favors a "standard" velocity 
law with of ca. 1-2 at most. 

5.2. Measured sub-peak acceleration from synthetic 
Ipv: O-star case 

In the previous section, we chose an inner radius for the 
LER compatible with the generally large emission line pro- 
file widths of WR stars. By contrast, in O-stars, emission 
line formation starts right at the base of the outflow, re- 
sulting in much narrower emission profiles, usually quite 
pointed rather than flat-topped. Thus because the whole 
wind is optically thin, a given emission line carries infor- 
mation on all wind heights, allowing us to constrain the 
inner wind acceleration from single-exposure line profiles 
(Sect. 3.2; Puis et al. 1996; Crowther et al. 2002), as well 
as from Ipv (Eversberg, Lepine, & Moffat 1998). 

To model the Ipv of O-star winds, we use the same ap- 
proach as in Sect. 5.1, this time setting the inner edge of 
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Fig. 11. Montage of synthetic Ipv based on radiation- hydrodynamics simulations of the line-driving instability for four 
optically thin emission lines whose LER extend out to 20i?» but start at different depths: 2.5 and 5.0 i?* (top left and 
right panels) and 10.0 and 15.0 i?* (bottom left and right panels). To enhance the visibility of migrating sub-peaks in 
the grayscales, the local maximum flux deviation is limited to 1%. 



the LER to very small wind heights, covering the values 
1.2, 1.5 and 2.0 R* (left to right panels in Fig. 11). The 
migration magnitude of profile sub-peaks is quite com- 
parable between all three simulations, and only slightly 
higher that for the situation where the LER starts at 2.5 
R* (Sect. 5.1). In our model (Fig. 10), although the wind 
acceleration magnitude increases towards the wind base, 



the wind ceases to be clumped below ca. 2 R*. Hence, no 
profile variations arise from wind emission at these depths, 
invalidating the use of our method based on Ipv to infer 
the wind acceleration. 

But, although the Ipv pattern is very similar in these 
three cases, the time-averaged line profile changes dramat- 
ically as the LER inner radius decreases. Indeed, the line- 
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Fig. 12. Same as previous figure but for 3 emission line profiles whose LER starts at 1.2, 1.5 and 2.0 i?* (left to right) 
but extends out to the same outer radius, i.e. 10.0 i?*. 



wing shape presents a very clear kink, at a profile location 
which corresponds to the wind velocity at which clump- 
ing appears. The emissivity of the Ipv diagnostics being 
density-squared dependent, it is suddenly enhanced when 
going from the unclumped inner wind regions through the 
clumped outer wind ones. Thus, if O-star wind clumping 
only started at a sizeable fraction of the terminal velocity, 
this feature would be observable. 



However, to our knowledge, such a profile kink has not 
been seen in O-star spectra which might be supporting 
the idea that their winds are clumped even in the inner 
wind region. The observed Ipv for the O supergiant star £ 
Puppis (Eversberg, Lepine & Moffat 1998) also suggested 
that the pattern of migrating subpeaks was compatible 
with a wind clumping starting right at the wind base. 
Moreover, Prinja ct al. (2001) reported a B/R asymmetry 
in the Ipv of HD 152408 (08:Iafpe or WN9ha), with the 
migration magnitude of profile sub-peaks twice as small 
in the red part than in the blue part of the line profile. 
In Sect. 3.2, we showed that disk occultation can lead to 
a significant flux deficit in the red part of the profile. But 
if this deficit is combined to a B/R dichotomy in Ipv, it 
also implies that the occulted innermost wind regions are 
also clumped. Thus, these various elements support the 
spectroscopic analysis of O supergiant stars (Hillier et al. 
2002, 2003; Crowther et al. 2002) which require that the 
wind clumping starts right from the base. 



6. Summary and future work 

We have reexamined the LER characteristics of WR Ipv 
diagnostics (Robert 1992). We find that, in general, these 
lines arc optically thin, showing a profile broadening that 
results from the finite extent of the LER rather than 
from turbulence. Furthermore, we find that they form 
much closer to terminal velocity than proposed in previous 
studies (LM), leading to a downward revision of the WR 
wind acceleration length scales f3R* from 20 — 80 Rq to 
10 — 20 Rq. Based on ID radiation hydrodynamics simula- 
tions of the radiative instability and the "Patch" method, 
we compute synthetic Ipv for Ciii5696A for WR111. We 
are able to reproduce well the observed Ipv for that star, 
although our model uses an O-star core radius. Both evo- 
lutionary and spectroscopic models predict small WR core 
radii, a few times the solar value, which thus makes WR 
wind acceleration length scales still considerably extended, 
although not to the extent proposed by LM. Forthcoming 
VLTI-AMBER observations will hopefully give an accu- 
rate evaluation of the spatial scales of WR stars and their 
associated winds, thereby setting the magnitude of the dis- 
crepancy between the inferred WR wind acceleration from 
observations and that expected from radiation driven wind 
theory. 

In our synthetic Ipv simulations for optically thin O- 
star winds, we find a clear signature for the onset velocity 
of wind clumping, which results from the strong sensitivity 
of recombination line emissivity to variations in wind den- 
sity. This introduces a kink in the line-wing profile shape, 
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Fig. 13. Optical emission lines of Ovi5290A, Ov5595A and Civ5808A observed in the Br93 spectrum. Superposed 
on each profile is a satisfactory fit with the model described in Sect. 3.3. The velocity centroid v e of the LER for the 
three modeled lines, going from left to right, are 40%, 65% and 85% of v^, while the velocity range of the LER is 
20% of Woo for each. No line opacity needed to be accounted for in order to reproduce the profile. Some turbulence 
is included: v tur b ~ 400 km s _1 , which corresponds to 10% of the terminal velocity. This illustrates why WO stars 
represent ideal candidates to reveal differential wind acceleration from Ipv. 



making the clumped outer wind regions relatively more 
emissive than the smooth inner wind ones. The absence 
of such a kink in even high-resolution spectroscopic ob- 
servations of O supergiants together with results from Ipv 
studies of such stars (Eversberg, Lepine, & Moffat 1998; 
Prinja ct al. 2001) suggest that indeed, wind clumping 
may well start right from the hydrostatic base in O super- 
giant stars. 

Finally, the WR- Ipv dataset obtained so far uses di- 
agnostic lines that form close to terminal velocity; such 
lines contain only poor information on the overall wind 
acceleration. Using optically thinner winds would yield a 
more valuable information on the wind acceleration. This 
is the case for OB objects but their wind density is how- 
ever so low that no line emission arises from above ca. 2-3 
R» (Hillier et al. 2002), providing no information on the 
outer wind dynamics. On the contrary, WO star winds 
have not only a low optical depth in the visual band but 
they also show strong optical emission lines of markedly 
different velocity widths, resulting from a strong wind ion- 
ization stratification (Crowther et al. 2000). Indeed, using 
the model of Sect. 3, we find that optical emission lines of 
oxygen and carbon form in regions located between 40% 
to 100% of the terminal velocity (Fig. 13). Thus, WO stars 
constitute excellent targets for Ipv studies to provide fur- 
ther observational constraints on the structure and dy- 
namics of WR star winds. 
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